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Abstract. We report results of a study of (integer) quantum Hall transitions in a single or multiple 
Landau levels for non-interacting electrons in disordered two-dimensional systems, obtained by pro- 
jecting a tight-binding Hamiltonian to corresponding magnetic subbands. In finite-size systems, we 
find that mesoscopic effects often dominate, leading to apparent non-universal scaling behaviour in 
higher Landau levels. This is because localization length, which grows exponentially with Landau 
level index, exceeds the system sizes amenable to numerical study at present. When band mixing 
between multiple Landau levels is present, mesoscopic effects cause a crossover from a sequence of 
quantum Hall transitions for weak disorder to classical behaviour for strong disorder. This behaviour 
may be of relevance to experimentally observed transitions between quantum Hall states and the 
insulating phase at low magnetic fields. 
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1. Introduction 

One of the central ideas used to analyse phase transitions between different integer quan- 
tum Hall states in two-dimensional electron systems is finite size scaling. Transitions be- 
tween plateaus in the Hall conductivity are characterized by a localization length £ which 
diverges at the critical point s c , £(s) ~ \s — s c \~ u . s can be energy E (in numerical 
calculations), or magnetic field B (in experiments). Finite size scaling asserts that vari- 
ous quantities for a finite system of size L will depend on the ratio L/£. (In experiments, 
the relevant size is the dephasing length, which is determined by temperature T through 
inelastic scattering mechanisms). Thus, a system realizes its finite size when L/£ < 1, 
when singularities associated with the phase transition are cut off by the finite size. This 
defines the range where finite-size effects become important. Therefore, long length scales 
always pose a challenge for numerical calculations, in which the system size is limited by 
the current computing power. Unfortunately, the integer quantum Hall transitions in higher 
Landau levels are among such cases. 
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To a large extent, the difficulties of identifying the nature of the quantum Hall transitions 
in weak magnetic fields come from the same origin. In this case, the mixing of neighbour- 
ing Landau levels exacerbates the situation even further. According to the global phase 
diagram, developed by early theoretical efforts [], for integer quantum Hall transitions, the 
Hall conductance can only change by e 2 /h at the transitions, whether in high magnetic 
fields or in low fields. Early experiments [] were consistent with this expectation (counting 
spin degeneracy). However, more recent experiments by Song et al. [], Hilke et al. [], and 
Kravchenko et al. [] cast doubts on the global phase diagram, as their data were interpreted 
as showing direct transitions from an integer quantum Hall state with Hall conductance 
ne 2 /h (n > 1) to an insulator. Subsequently, some numerical calculations [] appeared to 
corroborate the existence of direct transitions between quantum Hall liquids with Hall con- 
ductance of ne 2 /h (n > 1) to the low field insulator. However, it has never become clear 
whether such "transitions" are true quantum phase transitions or merely the manifestation 
of consecutive transitions which are too close to resolve. 

In this paper, we discuss the mesoscopic effects in the quantum Hall regime, in particular 
the effects of long length scales on the scaling behaviour and on the low-field transitions, 
using a tight-binding lattice model of two-dimensional non-interacting electrons in a ran- 
dom potential and a perpendicular magnetic field. We describe our model with truncated 
Hilbert space and various methods in Sec. ||] and ^respectively. Section ^] contains the 
discussions on the difference of the scaling behaviour in Landau levels with increasing in- 
dex, as a result of the increasing length scales. We then turn to transitions at low magnetic 
fields in Sec. WL before summarizing our results and conclusions in Sec. 



We consider a tight-binding Hamiltonian of non-interacting electrons on a two-dimensional 
square lattice with nearest neighbour hopping, subject to a uniform perpendicular magnetic 
field B: 



€i are the on-site random potentials, and the magnitude of the hopping matrix element t = 1 
is chosen as the unit of energy for convenience. As a consequence of the perpendicular 
magnetic field B, there is a magnetic flux <\> per unit cell given by 



where cf>o — hc/eis the magnetic flux quantum, and ao is the lattice constant. The random 
potential e^, are taken to be independent random variables, chosen from an identical dis- 
tribution (so ensemble averaged quantities are translationally invariant). We consider two 
specific distributions: 

(i) a rectangular box (i.e. uniform) distribution on [— W, W], 
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Figure 1. Density of states p(E) of a 14 x 14 square lattice with magnetic flux cpo/7 
per unit cell and W = 0.5 rectangular box distributed disorder. 

(ii) a Gaussian distribution with standard deviation a, 

p(e) = vfc ex >'(-^)' ,4) 

The calculation of the energy spectrum can be done for rational <fi/(f)o = p/q, where p 
and q are integers, chosen to have no common factors. The whole spectrum [] is commonly 
known as the Hofstadter butterfly, so named because of the butterfly-like self-similar pat- 
tern in terms of the parameter <p. For each pair of p and q, in the absence of disorder, the 
original tight-binding band is split into q magnetic subbands; these are then broadened by 
disorder. Figure |l] shows the density of states of a 14 by 14 square lattice with magnetic 
flux <po/7 per unit cell. The seven magnetic subbands are well separated by energy gaps at 
disorder strength W — 0.5. 

The Hall conductance of each of the q magnetic subbands, when they are separated by 
energy gaps, has been shown to be quantized []. The Hall conductance a r (in units of 
e 2 /h), when the lowest r subbands are occupied, is related to p and q by a Diophantine 
equation [], r = qs r + pa r , where s r is an integer chosen such that |a>| < q/2. The 
quantized Hall conductance has a topological origin, which is the Chern number of a fiber 
bundle defined by the magnetic Bloch wavefunctions on the torus, the magnetic Brillouin 
zone. For the case of p = 1 and q > 1 an odd integer, each of the (q — 1) side subbands 
carries Hall conductance e 2 /h. The center subband, on the other hand, has Hall conduc- 
tance of — (q— 1 ) e 2 /h, guaranteeing that the total Hall conductance of all the subbands of 
the tight-binding system is zero. The continuum model is achieved in the limit q — > oo, 
where each of the side subbands corresponds to a Landau level, while the center subband, 
with no counterpart in the continuum model, is mapped to infinite energy. For finite q, it is 
expected that the side subbands near the edges of the spectrum, correspond closely to the 
Landau levels in the continuum model. 

Just like in the continuum model [], we can project the Hamiltonian to a subspace 
spanned by the eigenstates in the subbands of interest, obtaining from diagonalizing the 
tight-binding Hamiltonian in the absence of disorder. This truncation process is equivalent 
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to the addition of a pseudopotential that pushes the subbands not of interest to positive 
or negative infinite energy. The projection can be easily justified in the strong field limit 
(weak disorder), where there are energy gaps separating the subbands. On the other hand, 
in the case of weak magnetic fields (strong disorder), numerical calculations show that the 
center subband widens with increasing disorder [] and the negative Hall conductance car- 
ried by the central subband neutralizes the Hall conductance of the opposite sign carried by 
side subbands. This has been interpreted as quantum phase transitions from the quantum 
Hall liquids (with arbitrary integer filling factors) to the insulating state []. However, since 
there is no center subband in the continuum model, we have chosen to work with a pro- 
jected Hamiltonian in a truncated Hilbert space with multiple side subbands only. This has 
a major advantage that unphysical effects due to the center subband in the discrete lattice 
model are avoided. (We note that the actual experimental system is on a discrete lattice, 
but that is on an atomic scale, while the physics we are concerned with is on the scale of 
the magnetic length, which is two orders of magnitude larger). We believe the projected 
Hamiltonian should keep the essential features of the low magnetic field physics, at least 
near the bottom of the spectrum where interactions with truncated higher Landau levels are 
weak. 



3. Chern Number and Thouless Number 



We consider a system of size L x x L y with generalized periodic boundary condition 

t(L x , y )\m) = e i6 *-y\m) (5) 

where t(L x y ) is the magnetic translation operator in x or y direction. The Hall conduc- 
tance of an eigenstate can be calculated by the Kubo formula as [] 

rr < m -ft () ■\- ie2fl V" ( m \vx\n}(n\v y \m) -h.c. 

a xy (m M ) - — ^ - . (6) 

Here \m), \n) are eigenstates for the boundary condition Eq. (^|) and v x , v y are components 
of the velocity operator. Thouless et al. showed that the averaged Hall conductance can be 
written as a quantized integral [] 



1 



T xy (m) = / d9 x d9 y a xy (m;9 x ,9 y ) = C(m) — , (7) 



where C(m) is an integer, known as the Chern number of the eigenstate. States with non- 
zero Chern number conduct Hall current and can be identified as conducting states []. The 
total Hall conductance a xy (Ef) for a given position of the Fermi level Ef is obtained by 
summing over the Hall conductances of the filled states. 

The longitudinal conductivity a xx , related to Thouless number, may also be used to dis- 
tinguish extended states from localized states []. Thouless number, which describes the 
sensitivity of an eigenstate to the change in boundary conditions, is defined in dimension- 
less form [] by the equation: 

9l(E) = ^ ~ - 2 a xx , (8) 
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where (SE) is the average shift in the energy level due to a change of boundary condition 
from periodic to anti -periodic, and (AE) = l/L 2 pi J (E) the mean energy level separation. 
In second-order perturbation theory, the change in eigenenergy is proportional to the longi- 
tudinal conductance []; the proportionality constant is of the order of unity, with a precise 
value that depends on the exact definition of the energy shift (SE) . 

As the size of the system increases, the Thouless number in the localized regime de- 
creases, while at the critical point the Thouless number goes to a universal value. In the 
scaling regime, the dimensionless quantity depends on system size only through the ratio 

9L (E) = ~g(L/0=9(L 1 ^\E-E c \), (9) 
where the localization length £ diverges at the critical energy E c , in a power law fashion, 
£(E)~\E-E C \- V . (10) 
We define the width of a Thouless number curve gL(E) as 

1/2 



W(L) = 



1 



W DOS (L) 



!Zo9L(E)dE 



(11) 



W{L) is normalized by W^os{L), the width of the corresponding density-of-states pl(E), 



Wnos(L) 



Sr^E-EcYpM dE 
JZPL(E)dE 



1 1/2 



(12) 



which is expected to be only very weakly dependent on L. Using Eq. (|lO|), the width of a 
Thouless number curve can be shown to vary with system size as a power law: 



W(L) ~ L~ l ' v . 



(13) 



Alternatively, we may compute the area A(L) under gi{E), as A(L) = gL(E)dE. 
In the scaling regime, we have 



A(L) - L- 1 ^ 
with the same exponent v. 



(14) 



4. Quantum Hall Transitions in a Single Subband 

4.1 The lowest subband 

We start with the simplest case - truncation of the Hilbert space to the lowest subband, 
which we compare with the numerous results in the lowest Landau level available to date. 
The critical behaviour of the quantum Hall transition in the lowest Landau level has been 
studied numerically by various groups []. Numerical results in different models agree that 
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Figure 2. Thouless number g(E) of the lowest subband with W/t = 5 for 10 x 10, 
20 x 20, and 30 x 30 lattices. The inset shows the normalized width of g(E) as a 
function of size L, on a log-log scale. 



the localization length diverges with critical exponent v ~ 2.35 ± 0.1 in the vicinity of the 
critical energy E c at the band center. In particular, Huo and Bhatt found the localization 
length critical exponent v = 2.4 ± 0.1, obtained from the finite-size scaling of the density 
of extended states identified by a Chern number calculation []. 

If we choose <f>/4>o = 1/5, the tight-binding band splits into five magnetic subbands, 
which carry Hall conductance of 1, 1, -4, 1, 1, in units of e 2 /h, respectively. We separate 
the lowest subband carrying Hall conductance e 2 /h from the rest of the four subbands (car- 
rying total Hall conductance ~e 2 /h) by truncating the Hilbert space as described above. 
We studied systems with size L = 10, 15, 20, 25 and 30, with a rectangular box dis- 
tributed random potential W/t = 5. We diagonalized between 200 and 1,800 samples 
with different random potential configurations from the ensemble, for calculating average 
quantities, depending on size. Figure [| shows the Thouless number of the edge subband 
for sizes L = 10, 20 and 30. We clearly see the shrinkage of the width of the Thouless 
number as the system size increases. Meanwhile, the peak value of the Thouless number 
for different sizes remains fixed around go ~ 0.21, which helps us locate the critical energy 
E c — —2.9. The width of Thouless number of different sizes calculated by Eq. ( pi] ) (and 
normalized by the width of the density-of-states), is plotted in the inset on a log-log scale. 
Fitting the width to a straight line (a power law), we can extract the localization length 
exponent v = 2.4 ± 0.1 according to Eq. (|l3|), which agrees with the exponent found for 
the lowest Landau level in the continuum model. The finite-size scaling fit suggests that in 
the thermodynamic limit (L — ► oo) the Thouless number, or the longitudinal conductance 
a xx , vanishes everywhere except at one single energy E c , where the localization length 
diverges. This demonstrates that, under truncation, the edge subband has the same scaling 
behaviour as the lowest Landau level. 

In order to test the universality of the scaling behaviour, we studied systems with vari- 
ous values of magnetic flux per unit cell, including 0/</>o = 1/3, 1/5, 1/6, 1/7, and 1/9, 
choosing either rectangular box distributed or Gaussian distributed potential, with disorder 
strength W/t (or a/t) varying from 1 to oo (t — > 0). We calculated the Thouless number 
of the truncated lowest magnetic subband of these systems up to L = 66 and averaged 
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over as many as 100,000 random potential realizations (for L — 3). With these choices of 
system configurations, we found the localization length critical exponent v = 2.3 ± 0.2, 
obtained from the finite-size scaling of the width of the Thouless number. The peak value 
<7o of the Thouless number was found for all different cases to be go — 0.21 ± 0.02. This 
indicates that v and go are universal in the scaling regime, which includes system sizes as 
small as L = 10 in most cases. 



4.2 The second lowest subband 

It has been recognized that the localization length, 

E-E r 



AE 



(15) 



in higher Landau levels is much larger than that in the lowest Landau level []. In the 
above equation, AE is the width of the density of states p(E). £q, the E'-independent 
prefactor, is distinct for each Landau level, and is of order of the lattice constant in the 
lowest Landau level. The localization length in higher Landau levels grows rapidly with 
increasing Landau level number (exponentially as its square in simple approximations), 
raising great difficulties in numerical calculation of the critical exponent in higher Landau 
levels. This difficulty has been reflected, in particular from the observation of non-universal 
(larger-than-expected) exponent, in early works []. The universal scaling behaviour can, 
therefore, only be restored in much larger systems, unless the localization length can be 
reduced. 

A potential with a non-zero correlation length is known to be able to reduce the local- 
ization length substantially, to a numerically acceptable value. Calculations for the second 
lowest Landau level suggest that the value of the localization length exponent v depends 
on the correlation length I of the random potential []. Mieck showed that v decreases from 
6.2 for Z = to 2.3 for / = Al c , where l c = (Hc/eB) 1 / 2 is the magnetic length []. Huck- 
estein also obtained for a correlation length I = l c the same scaling behaviour as in the 
lowest Landau level []. We found very similar behaviour in the truncated second lowest 
magnetic subband, which, again, shows the correspondence between our lattice model with 
the truncated Hilbert space and the continuum model. 

In our lattice model, a correlated potential is generated as follows. First, generate an un- 
corrected random potential e° with the desired distribution. We then introduce correlation 
between random potential on different sites by defining 



(16) 



where dij is the distance between sites i and j. The correlation function between the two 
sites 

(eiEj) ~e~4/< 2 , (17) 

with correlation length I, It should be emphasized that our correlated potentials have Gaus- 
sian correlation; to get an exponential correlation, an appropriate modification of Eq. ( |i"oj ) 
is necessary. 
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Figure 3. The width of Thouless number g(E) of the second lowest subband as a 
function of size L on a log-log scale. The inset shows the localization length critical 
exponent v calculated from the slopes of the different lines in main figure for different 
correlation lengths /. 



In this section, we study the Thouless number of the (truncated) second lowest subband, 
of a seven-subband system (4>/<fio — 1/7), with correlated potential generated using a 
rectangular box distribution for e^. After projection we set t = (infinite disorder limit) 
to reduce finite size effects. For sizes L = 14 to 49, we diagonalized 150 - 1,800 samples, 
depending on L. We calculated the Thouless number g(E) of the subband with various 
correlation lengths. Figure [3] summarizes the width of the g(E) normalized by the width of 
the density of states. On the double logarithmic plot, the width versus size follows straight 
lines apparently consistent with the scaling hypothesis Eq. (|l3|). However, the localization 
length critical exponent v extracted from the slope of the curves is not universal. The inset 
of Fig. |3| shows that v decreases from 6.8 with I = to 1.4 with I = 3, and then increases 
slowly to 1.7 with I = 5. This curve is very similar to the one obtained by Guo [] in the 
study of the critical behaviour in the first Landau level in the continuum model. The curve 
crosses the expected universal exponent v — 2.3 for I ~ 1. 

Figure^ also shows that W(L) is minimized with I = 2, suggesting that the localization 
length is reduced with increasing correlation length I of the correlated potential, at least for 
I < 2 . However, with I = 2, v = 1.5 is much smaller than the expected v — 2.3. This 
suggests that the random correlated potential introduces an irrelevant length scale (, which 
depends on the correlation length /, so that 

9l{E)= 5(7=7)= (18) 

Only when L, £ ^S> ( can we restore the universal behaviour of quantum Hall transitions. 
Due to this additional length scale £, Thouless number, in general, does not have either the 
universal functional form or the universal peak value go in finite systems. 

We might expect the irrelevant length scale £ to be less important as the localization 
length £ increases and dominates when the band center (E c = 0) is approached. On the 
other hand, the irrelevant variable clearly affects our results in the tails where £ is small. 
The width of the Thouless number defined by Eq. (pT[), however, depends strongly on the 
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Figure 4. The localization length critical exponent v in the second lowest subband, 
calculated for different correlation length I of the random potential, from the scaling of 
the area A(L) with size, L. The inset shows the area A(L) of Thouless number g(E) 
of the second lowest subband as a function of size L on a log-log scale. 

scaling behaviour of the tail, since the value of g(E) in the tails is amplified by multiplying 
(E — E c ) 2 in the numerator. Therefore, the localization length exponent obtained from the 
finite-size scaling of the area A(L) of the Thouless number may be more reliable (which 
includes the finite-size effect of the peak conductance as well). 

Figure |] inset shows the area of Thouless number g(E) of the second lowest subband 
as a function of size L. On the log-log plot, the area A(L) versus L follows straight lines 
with different slopes. The numerical value of the localization length critical exponent v 
calculated based on Eq. ( |l4| ) is summarized in Fig. fij The universal value v = 2.2 ± 0.2 
is restored for I > 1. These values of I are consistent with the observation in Fig. || that 
W(L) is minimized with I = 2. v = 6.5 with i = is, again, due to the finite size effects 
as a result of the large localization length. 



4.3 The third lowest subband 

In even higher subbands (Landau levels), the localization length £ is expected to be even 
larger than £ in the second subband, with both uncorrelated and correlated potential. We 
confirmed, in our calculations, that the effects of Gaussian correlated potential on the scal- 
ing behaviour of the third lowest (side) subband of the seven-subband system agree quali- 
tatively with that of the second lowest subband. Once again, to minimize the effect of the 
irrelevant length scale £, we plot the area A(L) of the Thouless number curve versus the 
system size L in the inset of Fig. ||. For L > 21, A(L) follows straight lines on double 
logarithmic plot. Figure [| summarizes the localization length critical exponent v calcu- 
lated according to Eq. (|1J). v decreases from 16 with I = to v ~ 2.3 with I > 2. The v 
versus / curve is very similar to the curve in Fig. ^, except that, in the second subband, v 
is smaller with I = 0, and drops to v ~ 2.3 with 1 = 1. This is consistent with the fact that 
the localization length £ is smaller in lower subband. 

Therefore, due to rapidly increasing localization length in higher Landau levels and the 
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Figure 5. The localization length critical exponent v calculated at different correla- 
tion length /. The inset shows the area of Thouless number g(E) of the third lowest 
subband as a function of size L on a log-log scale. 



consequent mesoscopic effects, finite size scaling of the area A under the Thouless number 
curves seems to be more reliable than the width of the Thouless number curve. A universal 
exponent v ~ 2.3 can be observed for correlated potential with certain correlation length 
in finite size samples. It is worth emphasizing that the universality, including localization 
length exponent and critical conductance, should be observed in large enough samples. 
These are beyond our current computing capabilities, but should become more accessible 
as computers get more powerful. 



5. Quantum Hall Transitions in Multiple Subbands 

In this section, we briefly present our results for multiple subbands. More details can be 
found elsewhere []. In a lattice with <fi/<f>o = 1/7, we projected the Hamiltonian to the 
lowest three subbands, carrying a total Hall conductance of 3e 2 //i. We found that, with 
increasing disorder, Thouless number g(n), calculated as a function of filling factor n, 
evolves from three peaks separated by gaps at integer fillings to a smooth curve with no 
dips. Meanwhile, the plateaus of Hall conductance a xy at integer fillings, characteristic of 
quantum Hall transitions, disappears. 

Figure ^ shows g(n) and a xy (n) as a function of filling n, with disorder strength W/t 
= 4.5, which corresponds to the strong disorder (weak field) limit. We calculated g(n) in a 
28 x 28 lattice, and a xy (n) in a 14 x 14 lattice. The dip in g(n) near the upper end of the 
spectrum is affected by the truncation of higher subbands. However, at the lower end, we 
believe our data shows the correct low-field physics, since the energy levels in the subbands 
not included in the truncated model interact only weakly with the low-lying energy levels. 
For n < 1.5, we find that g(n) can be fit by a straight line, while the Hall conductance a xy 
is much smaller than the diagonal conductance, and is essentially zero. Such observations 
are generically available for systems of different magnetic flux and number of subbands. 

We attribute our results to the "classical" behaviour of the two-dimensional electron 
gas. The linear dependence of the longitudinal conductance on electron concentration is 
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Figure 6. Thouless number and Hall conductance (in units of e 2 /h) in the truncated, 
three lowest subbands with W/t = 4.5. 



consistent with the classical Drude conductivity a xx ~ n. This result seems to suggest that 
for finite sizes, due to the presence of long length scales, quantum interference effects can 
be cut off at strong disorder, or equivalently at weak magnetic field. This could explain the 
experimental phase diagram []. The concurrent vanishing a xy (lower panel of Figure ^ 
supports this scenario too. For further details, the reader is referred to Ref. []. 



6. Conclusions 



In summary, we have studied a tight-binding lattice model for non-interacting electrons, 
using projection to enable truncation of the Hilbert space to different combinations of 
the magnetic subbands. In a single magnetic (side) subband, we obtained, through 
finite-size scaling of Thouless number, the critical exponent v of the localization length 
£ ~ \E — E c \~ v , as v = 2.4 ± 0.1, in good agreement with previous numerical results in 
continuum lowest Landau level model. The localization length critical exponent is univer- 
sal, regardless of the strength and distribution of the random potential or magnetic field. 
In higher Landau levels, however, the rapidly increasing length scales can alter the scaling 
behaviour in finite size systems, a manifestation of mesoscopic effects in the quantum Hall 
regime. With the introduction of correlated potential, we demonstrated that the localiza- 
tion length can be reduced, leading to universal scaling behaviour of appropriately chosen 
quantities. 

In a truncated system involving multiple magnetic side subbands (corresponding to mul- 
tiple Landau levels), we found that the Hall conductance a xy of finite lattices evolves from 
a series of well-defined plateaus at ne 2 /h for various integer Landau level filling to a sin- 
gle sharp rise at the upper end of the energy spectrum. At the same time, the Thouless 
number of the (finite) system become linear in filling factor n. The linear dependence is 
interpreted as a crossover to the classical (unquantized) Hall response, as the system size 
becomes comparable to, or smaller than, the very (exponentially) large localization lengths 
of the electronic eigenstates in this limit. While the sizes studied here are considerably 
smaller than those corresponding to experimental samples, our findings raise the question 
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of the proper interpretation of the experimentally observed low field transition, and the 
possibility of the influence of a quantum-classical crossover in experiments. 
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